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Abstract: We derive the corrections to the Friedmann equation of order in the Randah- 
Sundrum (RS) model, where two 3-branes bound a shce of five-dimensional Anti-deSitter 
space. The effects of radion stabilization by the Goldberger-Wise mechanism are taken into 
account. Surprisingly, we find that an inflaton on either brane will experience no order 
corrections in the Hubble rate H due to its own energy density, although an observer 
on the opposite brane does see such a correction. Thus there is no enhancement of the 
slow-roll condition unless inflation is simultaneously driven by inflatons on both branes. 
Similarly, during radiation domination, the correction to -ff on a given brane vanish 
unless there is nonvanishing energy density on the opposite brane. During the electroweak 
phase transition the correction can be large, but is has the wrong sign for causing sphalerons 
to go out of thermal equilibrium, so it cannot help electroweak baryogenesis. We discuss 
the differences between our results and exact solutions in RS-II cosmology. 
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1. Introduction 

The idea of our world being a 3-brane embedded in a higher dimensional bulk has provided 
an attractive explanation for the weakness of gravity compared to the other fundamental 
forces ||^, ^ ^ . In this framework, 4-D gravity is weak at distances much greater than the 
fundamental Planck length, usually assumed to be the TeV~^ scale, but it quickly becomes 
strong at shorter distances. This idea has rich consequences for LHC physics, and it also 
implies modifications to cosmology at temperatures approaching 1 TeV. Although it may 
not be feasible to test deviations from standard cosmology which happened so early, this 
scale is sufficiently close to the electroweak phase transition to give some hope of observable 
consequences, for example to electroweak baryogenesis. Moreover the nature of inflation 
can be significantly altered by extra-dimensional effects [Q, ^, 18]. 



A seminal work on brane cosmology [Q] observed that the Priedmann equation on 
a 3-brane embedded in a flat 5-D spacetime takes the form = g^g^ instead of the 

standard form = Here M5 is the 5-D Planck mass and p is the 4-D energy 

density on the brane. This strange result arises from the fact that the brane is assumed to 
be infinitesimally thin compared to the size of the extra dimension (y), which introduces 
a discontinuity in the derivatives of metric elements with respect to y at the position of 
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the brane. The dependence is strongly contradicted by big bang nucleosythesis. It 
was subsequently observed ^ that by adding a negative cosmological constant A5 to 
the 5-D bulk, and a tension T to the brane, the Friedmann equation generalizes to = 
6M? + ^^^^^ ^^^^^ = ^^{l + ©(-^fr)) if one chooses = -6A5MI 



and = In other words, we recover the standard Friedmann equation, modified 



-■5 

18M| ~ ~3~ 

by corrections which are small as long as p ^ T ~ TeV'*. In addition, the two relations 
between T and M5 are equivalent to those needed for the Randall-Sundrum solutions, 
where the bulk is a region of 5-D Anti-DeSitter space. 

Although the problem of how to recover the correct Friedmann equation was partially 
resolved, refs. § pointed out that in the two-brane model of |^ (RS-I), could be real 
only if p was assumed to be negative on the second brane, whose tension is negative. In 
addition one had to assume a fine-tuned relation between the energy densities on the two 
branes. This untenable result was especially troubling because the negative tension, or TeV 
brane, is the one where we are assumed to be living in order to perceive 4-D gravity as being 
weak. Ref. [10| pointed out that these problems could be avoided if the Tss component of 
the bulk stress-energy tensor had the right properties. It was subsequently shown ||ll|, |l2| 
that stabilization of the size of the extra dimension (the radion mode) was necessary and 
sufficient for recovering acceptable cosmology on the TeV brane. The unphysical require- 
ment of a fine-tuned negative value for p was merely a consequence of assuming the extra 
dimension was static when generically, in the absence of any stabilizing force, it would tend 
to expand or contract. The positive effects of stabilization were subsequently confirmed 
in greater detail by references [13, 14] and [|l^ ]. The problem of losing stabilization at 
temperatures above the TeV scale and tunneling back to the stabilized ground state was 
explored in ref. |16|. 

The issue of radion stabilization exists only in the RS-I and ADD solutions; in 
RS-II 1^ the extra dimension is infinite in size, there is only a single brane, and there is 
no radion — its wave function diverges away from the brane, and so it is not a normalizable 
mode if there is no second brane to compactify the space. Because of this simplicity, much 
of the early work on brane cosmology has focused on the RS-II solution. In this case the 
above-quoted form of the modified Friedmann equation becomes exact, 

ii' = '-^f{i + S,) (1.1) 



3 V 2TJ 

as long as higher-derivative corrections to the gravitational action are neglected^ and the 
bulk is empty. Perhaps the most important application of the p/T correction is to increase 
the expansion rate over its normal value during inflation |^, . This allows for the possibility 
of fulfilling the slow-roll condition with steeper inflationary potentials than would normally 
be admissible, with the possible observable consequence of a spectral index of density 
perturbations which is smaller than in conventional models of inflation. It should be kept 
in mind that the effects of the p/T correction should not be trusted quantitatively whenever 



^ Since the discontinuity in the derivative of the warp factor at the brane is proportional to p as a 
consequence of the 5-D Einstein equations, one may expect that higher derivative terms will bring higher 
powers of p. 
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they become large, precisely because of the neglect of higher derivative operators which 
are only suppressed by the energy scale ~ T^/^. Beyond this scale one really needs the 
full quantum theory of gravity, presumably string theory, to make quantitative predictions. 
Thus the p/T corrections should be regarded as indicative of the kinds of new qualitative 
effects that could be expected at energy densities above the quantum gravity scale. 

More recently, the cosmology of the RS-I model has attracted renewed attention. The 
AdS / CFT correspondence has been used as a way of better quantifying the quantum gravity 
effects: Kaluza-Klein excitations of the bulk graviton modes are supposed to be equivalent 
to bound states of a strongly coupled, nearly conformal field theory residing on the TeV 
brane [^] . At temperatures above the TeV scale, the TeV brane is supposed to be hidden 



behind a horizon [18| associated with a black hole which formed in the bulk [|T^. It has 
been noted that the emergence of the TeV brane from the horizon may occur around the 
same time as the electroweak phase transition |2C]. Ref. [16| provided an alternative picture 
of a first order phase transition which leads to the appearance of the TeV brane at this 
epoch. 

There have been numerous investigations of the effects of a bulk scalar field on brane- 
world cosmology ||2l|. In the present work we shall be concerned with the 0{p'^) corrections 
to the Friedmann equation which arise from the 5-D nature of the underlying geometry 
in the RS-I model, using the Goldberger-Wise (GW) [^2| mechanism to stabilize the size 
of the extra dimension. In general one no longer expects eq. ( |1.1| ) to hold in the RS-I 
model because of the additional effect of the dynamical radion. The radius is displaced by 
the cosmological expansion, which changes the strength of gravity and thus back-reacts on 
the expansion rate. Although one should not trust these corrections quantitatively when 
they start to become large (for the reasons discussed above), they accurately predict the 
deviations from standard cosmology as one starts to approach the quantum gravity scale 
around 1 TeV. 

The plan of the paper is as follows. In section 2 we expand the 5-D Einstein equations 
to second order in a perturbation series in p and p*, the energy densities on the TeV and 
Planck branes. In section 3 the equations are solved to find the ingredients from which 
one can infer the effective 4-D Friedmann equations. In section 4 we present the results for 
the Hubble rate and the acceleration as measured on either brane, for arbitrary equations 
of state, and we discuss the implications for cosmology. We also compare the results to 
where there is only a single brane or where there is no bulk scalar. The reader who is not 
interested in the details of solving the Einstein equations can go directly to this section. 
Section 5 gives a summary and conclusions. Some technical details can be found in the 
appendices. 



2. O(p^) Einstein Equations 



Following closely the formalism of ref. [15|, we will look for cosmological solutions to 5-D 
gravity coupled to the Goldberger-Wise scalar field $, which is responsible for stabilizing 
the radion, and to two branes: the Planck brane, located at yo = 0, and the TeV brane at 
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Ui = 1. The action is 

+ j d^X^iCmfi - \y=0 + J d^X^iCm,! " ^l($)) \y=l, (2.1) 

where ac^ = M^^, and the potential V{^) is left unspecified for now. The brane contribu- 
tions are a sum of matter, represented by 



^m,0 ~ P*', 



-mA 



(2.2) 



and tension, which is the value of the brane's scalar potential Vi($(?/i)). The matter 
Lagrangians cannot be written explicitly for cosmological fluids, but their effect on the 
Einstein equations is specified through their stress-energy tensors ( |2.7| ). Our ansatz for the 
5-D metric has the form 



ds^ = n\t, y) df - a\t, y) ^ dxj - b\t, y) dy'^ 

i 



(2.3) 



We will make a perturbative expansion in the energy densities p, p* of the branes around 
the static solution, where p = = 0: 

N{t,y) = Ao{y) + 6Ni{t,y)+6N2{t,y); A{t,y) = Ao{y) + 6Ai{t,y) + 6A2{t,y) 
b{t,y) =bo + 6bi{t,y)+6b2{t,y); ^{t,y) = Mv) + S^i{t,y) + 6^2{t,y). 

(2.4) 

The subscripts on the perturbations indicate their order in powers of p or , both of which 
are taken to be formally of the same order for the purposes of developing the perturbation 
expansion, although in actual order of magnitude we may consider them to have very 
different values. This ansatz is to be substituted into the scalar field equation, 

dt (^^a^^) - dy Qa^n$'^ + a^n [bV + Vo5{y) + Vl5{y - 1)] = 0, (2.5) 

and into the Einstein equations, Gmn = K-'^Tmn- Here and in the following, primes on 
functions of y denote while primes on potentials of $ will mean The components 
of the Einstein tensor which do not vanish identically are 



G, 



00 



a\ ^db "-^ [ _|_ I 



a 



+ — 



a J ab b"^ 
l\ 2 



+ 2- 



a n b n 



a n 



a n 



n 



a'b' 
ab 



b' a' a" n" 
2-- + 2— + — 

b a a n 



- +2 2- + - -2- + - 



a n 



a n 



a n 



-4- 



Go5 

G55 



n' 


a 


+ 


a' 


b 










n 


a 




a 


b ' 












a' 


( 


a' 


+ 


n' 


a 




a 


n 



a 
a 



62 



a a 



a \ a 



n 



n 



(2.6) 



and the stress energy tensor is T„ 
On the branes, is given by 



gmn{V{^)+k)+dm^dn^ - ^d^m^9mn m the bulk. 



+ 



b{t,0) 

S{y - 1) 
b{t,i) 



diag(Vb + P*, Vq - p^,Vo - p^,Vo - p*, 0) 



diag(Vi + p,Vi- p, Vi -p,Vi- p, 0) 



(2.7) 



(Later we will assume that the potentials Vq and Vi are very stiff and are vanishing at their 
minima, so they can be neglected.) 

The perturbative solution has been given to first order in p, p^, in [^]; here we want to 
extend the analysis to second order. Following |15|, it is useful to work with the following 
linear combinations of the metric tensor perturbations: 

*2 = SA'^ - A'o^-^ - y $;^<5$2 - ^ (^5^[ + <^'J-h}j s<^,- T2 = dN^ - SA'^. (2.8) 

These are convenient variables because they appear in a natural way in the boundary 
conditions at the two branes. By integrating the (00) and (ii) Einstein equations in the 
vicinity of either brane, we obtain the perturbed Israel junction conditions, 



^2{t,0) = + -rbop,-^ 
6 60 



T2(t,0) 



t;y=0 

k2 5bi 



6 60 



t-y=l 



t;y=0 



T2{t,l)= + —bo{p+p) 



(2.9) 

(2.10) 



t;y=l 



The analogous quantities at first order were found to be 

= SA'-A'J^- Ti = 5N[ - 5A'^ 

Oo 6 



(2.11) 



in [15 1, and they satisfy the same boundary conditions as in ( 2^ , 2.1C| ), but with the re- 
placement ^ — > 1. 

The Einstein equations take a form resembling ordinary second order differential equa- 
tions in y for the variables ^2 and T2, where all the time dependence of the latter is im- 
plicit, through /9(t), p*(t), p{t) and p^.{t). The quantities ^ and ^ function like constants 
of integration in this context, whose values are fixed by imposing the boundary conditions 
(2.E, 2.1C| ). This procedure results in o.d.e.'s in time for ao{t), which are the Friedmann 
equations we are seeking. Accordingly, we need to expand | and | in powers of p. Recall 
that in ordinary cosmology, (^)2 ~ ^ ~ p. Therefore each time derivative of oq counts as 
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half a power of p. In the brane cosmology we have (^)^ ~ ^ ~ P* + 0{pl, pp^, p^), so we 
expect that 



Oo 
2 



SvrG 



[p^ + n^p] l + 0{p^,p) 



^1 +1- 



ao 



ao 



+ ... 



00/(1) V«0/(2) 



(2.12) 



where we have used subscripts to denote the power of p or p^, contained in the associated 
term. Here = e""^"*-^-* is the value of the warp factor at the TeV brane, which solves the 
hierarchy problem in the RS scenario, and the linear combination p* + Q^p is the result 
that has been proven [11, ^ to appear in the lowest order Priedmann equation when the 
radion is stabilized. We therefore expand ^ and | in the following manner: 




+ o{p'y, 



(2.13) 



(2.14) 



In terms of these variables, we can write the second order Einstein equations (specifi 

2 

cally, the linear combinations (00), (00) + ^{H), (05) and (55), in that order) as 



ao 



(2) 



2 r:^ 



ao 7(2) 



O0/(2)/ 



Oq 



Ts + ^-2 + To5 



(2.15) 

(2.16) 
(2.17) 



«0 7(2) 



"^0 7(2)^ 



-) +(-) l^ge2^"=^[,(4^2 + T2) + y($'o'<5$2-$[,<^$'2 + 'J>^?^)+-^55 



(2.18) 



and the scalar field equation as 
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(2.19) 



where all the dependence on first order quantities squared is contained in the functions 
J^\i,, J^r, J^05, -^55 and which can be found in appendix A. 

We have used the symbolic manipulation features of Matlab and Maple to arrive at this 
and the other results in this paper, giving us confidence that there are no algebraic errors, 
despite the complexity of the formulas. We emphasize that all the terms with subscripts 
"1" (or ^) are already explicitly known from the solution given by ref. [^]. The unknown 
quantities, with subscripts "2" (or |) thus constitute a relatively manageable part of these 
complicated looking equations. Having thus derived the second order equations, we are 
now ready to find their solutions. 



3. Solutions 



We wish to solve the perturbed Einstein equations ( ^.15 - 2.18 ) for the unknown quantities 
(do/ao)^2) and (ao/ao)(2)) which are needed to find the corrections to the Friedmann 
equations on the TeV brane. To simplify this task, we will henceforth work in the limit of 
stiff brane potentials, Aj oo, which ensures that = at each brane. This allows us 
to choose a gauge in which vanishes everywhere in the bulk and to solve explicitly 
for 



^ = ^ [^'(P - + (p. - 3p*)(G - A',e^^^')] 



where 



G{y) 



-2Ao 



dy 



lg2Ao(s/) _^ ^^e4Ao(s/) 
The other terms which are first order in p were also found in 



~2Ao 



dy 



to be: 



^-1 
Ti 

where 

and 



6(1 -f]2) 
2(1 



„4Ao(j/) 



e^Aoiy) (_ir(y)(5^4(^ + p)+p,+ p.) + {n\p +p) + n\p, + p,))) 



F{y) = l-{l-n 



Ao(l) 



(3.1) 
(3.2) 

(3.3) 

) 

(3.4) 
(3.5) 
(3.6) 



is the warp factor evaluated at the TeV brane. 

^In fact, the procedure of using a small coordinate tranformation to set = can be reiterated to 
set 5^i{t,y) = to all orders, so that all effects of p and p* on the bulk stress-energy component T55 are 
contained in bi(t,y). 
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We recall that the zeroth order static solutions are approximately (ignoring higher 
order in Vq corrections) | 



My) = voe-^^'^y- 

where we have normalized ^o(O) = 0, and introduced 



Ao{y)^kboy + -vl{e-''^''''^y 



4 + 



fc2 



m 

W 



Here 



(3.7) 

(3.8) 
(3.9) 



k = V-k2A/6 

is the inverse AdS curvature scale, which should be somewhat below the Planck scale so 
that higher derivative operators like do not invalidate the solution. 

The above solution is an approximate one in the limit of small vq for the model con- 
sidered by GW [22 1, where the bulk scalar $ has only a mass term. It was shown by ref. 
1 23] that this solution is actually exact when $ has a quartic interaction which is related 
to the mass in a certain way. We are therefore justified in considering ( |3.7| ) to be exact, 
for the appropriate choice of V{^), in all that will follow. If on the other hand we decided 
to use the pure m?^'^ potential for V{^), the solutions ( |3.7| ) receive corrections which are 
a power series in K^foe"^'^'^'"'^: 



Ao{y) 
My) 



kboy + 



12 



(^~2khoy 



1) + 



96(1 + e) 



-Akboy 



M _^~3kboy _^ 



(3.10) 



12(1 + e) 

Here (pQ differs from vq by terms of order k^Vq, to satisfy the boundary condition (/)o(0) = vq. 
We shall find that none of our main results depend on which bulk scalar potential is used. 

As was previously shown, the first two Einstein equations yield differential equations 
for ^2 and T2. They are solved by: 

2 



^2 


= g4Ao{?/) 




y^e-4Ao(s/) 1 










T2 


= g4Ao(j/) 


Cr + 














dy 



(3.11) 



^ \ble^^^^^y)-Fr{y,t)\dy 



(2), 



(3.12) 



We do not require the $ equation of motion, since this is derivable from the Einstein equa- 
tions. The boundary conditions ( p.S|j2.ld ) allow us to eliminate the constants of integration 



C41 and Cx to find (ao/ao)(2) and (do/ao)(2)- 
2 



,«oy(2) 



«0 7(2) 



1 



bUie-^^^dy 
1 



«oy(2) 2blJ^e-^^ody 



^-2(0) - 0*^'2(1) + 

0^X2(1) -T2(0)+ / 
Jo 



-^^'^'^y^T^dy 
,-^My)jr^dy 



(3.13) 
(3.14) 
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Notice that the terms ^'2(0), ^'2(1); '^2(0), '""2(1)) are given in terms of p, p*, p and by 
the jump conditions ( |2.9[j2.1C| ). We will denote the equations of state for matter on the 
branes by 

p = ujp, p*=a;*p*. (3.15) 
To evaluate the integrals of J-^ and J-t in ( 3.13| - 3lT^ ), we need the first order solutions 



"fi, Ti, given by (^-3^), as well as 6A1 and 5Ni. These can be obtained from the 
definitions ( |2.1lD and the choice of gauge 6^1 = 0, as dAi{y) = Jq{Aq^ + ^i{y))dy and 



6Ni{y) = 6A1 + Jq Ti{y) dy. Moreover, we need (ao/«o)(i/2)) given by ( 2.12 ), and time 



derivatives of the first order perturbations, 6A1, SNi, 5bi. To reexpress the latter in terms 
of p and p*, we use the following relations, which can be derived from the (05) Einstein 
equation 

p = -3^(1 + u)p + 0{p^/^) ^ -v'247rG(p, +17V)(1 + ^) P 



p, = -2,Al + oj^)p* + 0{pl'^) ^ -V247rG(p* + rJV)(l + w*)p* 
p ^ 12^G(1 + uj)p [3(1 + uj)n^p + (3 + + 2l<;)p*] 

p, ^ 12^G(1 + w*)p, [3(1 + w*)p* + {3 + uj + 2a;*)pO^] (3.16) 

where we have assumed that uj and are constant in time. (This assumption would only be 
important during a period of transition such as going from radiation to matter domination.) 
In these expressions, the 4-D effective Newton's constant is defined by integrating the 5-D 
gravitational action over the extra dimension to obtain 

SttG = {2hQ e-'^^Hi^ . (3.17) 

We now have general expressions allowing us to obtain the second order corrections 
to the Friedmann equations. However the integrals appearing in these expressions (of the 



form j gConst.xAo^y^ cannot be performed analytically for ^0 given by (3/7). To overcome 
this difficulty, we will need to make a further approximation, namely that the radion mass 
is small. In the small- limit, the radion mass was found to be 

ml - -^2^,2,2^2^^2+2. (3_-Lg) 

3 



in ref. |14], which also agrees up to factors of order unity with the result found by ref. 
1 24]. By expanding in the small parameter we can express the warp factor 

as e^^o = e^(^''os/+'?[e-2''=''o^-i]) = ^ckh^y ^(c^fg-a^fc^o?/ _ i])«. The integrals can then be 
performed exactly, order by order in powers of v\. In fact, the appearance of I/^^q^ in ( |3.1| ) 
means that the leading order term in this expansion will actually be Oiv'^'^y). 

At this point we are ready to combine the terms like (do/ao)^2) ^"^^ (^o/fio)(2) with 
those coming from (ao/ao)(i/2)'^^i and bA\^ to obtain the complete results for d/a and 
d/a, as dictated by eqs. ( |2.13| , 2T^) . However, these are not yet the physical Hubble rate 



nor acceleration when evaluated on the TeV brane, due to a further correction from blSi\ 
that must be applied, and which will be explained in the next section. Because bA\ and 
(5iVi depend on y, the rates a/ a as measured by observers on the two branes will differ. 
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4. Physical Friedmann Equations on the Branes 



Although the main results found above, ( 3.13D -( p.l4 ), will resemble the Friedmann equa- 



tions when expressed as functions of the energy densities on the branes (after being com- 
bined with the appropriate corrections from 6Ai and 5Ai), they are not yet written in 
terms of the standard Friedmann-Robertson- Walker (FRW) time variable for TeV brane 
observers, which we shall denote by r. On the TeV brane, the lapse function 'n?{t, 1) has 
also received corrections of the kind we are interested in. Therefore we should transform 
to r using 

dr = ^^M^ dt = ^-SN,{t,i)-SN,it,iy- ^4^1) 
no{t, 1) 

where no{t, 1) = ao{t, 1) = 17 is the warp factor without any perturbations due to matter. 

At first it may not be obvious why we should use only the perturbations, 5Ni, to relate 
r to t instead of using the full warp factor. Of course, one is always free to rescale the time 
coordinate by a constant factor; this cannot change any physical observables. There are 
several ways of seeing why the above choice is the most straightforward one. First, it gives 
the same time coordinate as that in which the resolution of the gauge hierarchy problem 
was couched in the original RS paper [Q]. Specifically, using the time coordinate as defined 
above, in the absence of the cosmological perturbation, the action for a free scalar field on 
the TeV brane has the form 

= i|<i^(((vi^)-(Vvpn-17W) (4.2) 

where (p = This was the argument used to show that a bare mass of uiq ~ Mp 

translates into a physical mass of m = ^IrriQ. If on the other hand we had defined dr = Qdt 
as we might have been tempted to do in ( |4.1| ), we would get 

S ^\ j drd^xn^ ((02) - J7-2(V0)2 - mg02) 

= dTd\{{ip'^)-n-^{V^f -mlif'^) (4.3) 

where now (p = ^l^^'^cp. In these variables it is no longer obvious that physical masses are 
suppressed relative to the Planck scale. The resolution is that one must also reconsider 
the 4-D effective gravitational action in the new time coordinate. The result is that the 
4-D Planck scale is exponentially larger than the underlying 5-D gravity scale, (SvrG)"^ = 
^l^^K^"^ /k. In this picture, the fundamental scales are all taken to be of order TeV, while 
the 4-D Planck scale is enhanced by Q^^. Another, simpler, way of deriving the same 
result is to do everything with the warp factor defined to be unity on the TeV brane, and 
large on the Planck brane. In either method, the ratio of physical particle masses to the 
4-D Planck mass is the same. 

The consequence of all this is simple: to find the physical Hubble rate at some position 
in the bulk y in terms of the complete expression for d/a, one merely multiplies by dt/dr = 
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JNi+SN2\ . 

e \y. 

H^{l + 5Ni)- (4.4) 
a 

Since the leading contribution to H is already 0(p^/^), we can ignore the 6N2 corrections; 
these would be 0{p^/'^). Similarly, the second Friedmann equation, which goes like if, 
becomes 

ar at 

/ ■■ / ■ \ 2\ 



^(l + 2J7Vi) -- - \+6Ni^ (4.5) 
ya \aj J 00 

to the order we are calculating. All dots are still time derivatives with respect to t, not r. 

The general expressions for the order corrections to and dH/dr are rather 
cumbersome if expressed at an arbitrary position in the bulk. Here we give the results for 
H evaluated on either the TeV or the Planck brane, to leading order in the radion mass 



squared, mf, eq. (|3.18| ). We have defined p = p since this is the physically observable 



energy density on the TeV brane, whereas p is the bare value. For each term pi, p^^p, fP, 
we give only the leading dependence on the warp factor 17, which is usually assumed to be 
<C 1 to solve the hierarchy problem. Thus we have found the leading corrections to the 
Hubble rate and acceleration in a simultaneous expansion in p, p^, l/rrij. and Q,?' 

n, SvrG / 27rG' / , , , ^ n 

+4(1 - 3a;,)(4 + 2,u^)^^pl + 4(1 - 3l^)(4 + 3c<j)pp,))^ (4.6) 
H\=, = ^-^[-p + p^ + (9(1 - 3u;,)(l + u:,)pln^ - (1 - 3^)(7 + 



dH 
d7 



+2Q2p/),[2(l - 3w)(2 + 3uj) - 3(1 - 3w*)(l + w)] ) j 
-AttG ( p(l + u) + p^{l + w*) 



(4.7) 



dH 
d7 



y=l 

+ + oj,){l - 3w,)(13 + 9uj,)pI + 9(1 - 3cu)(l + coff 

- (1 - 3w)(2(l + w,) + 2(1 + w) + 6(1 + Lof + 3(1 + w)(l + uj^))pp*] 

= -AttG ( p{l + u) + p^{l + w,) 
y=o V 

-I- 



(4.8) 



+ ( -4(1 + ^)(1 - 3^)p2 + + ^^)2(i _ 3^^)f^4^2 



^We have factored these expressions to make it more clear that they vanish when lo = lu, — 1/3. The 
asymmetrical appearance of terms like 2(1 — 3lj)(2 + 3uj) — 3(1 — 3lj,)(1 + uj) is not due to a typographical 
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+ [6(1 - 3uj){l + ujf + 2(1 - 3uj){l +uj)- 2(1 - 3uj){l + w*) 

- 4:{l-3iu^){l + uj)+3{l + u,)il-3u){l+u)]p^p]] (4.9) 



These, finally, can be considered to be the main results of this paper, except in the case 
of u; = = 1/3. As can be seen, in that case the quadratic corrections vanish at order 
1/m^, so we have to go to higher order in to get the leading result for radiation, as will 
be described below. Two comments are in order. (1) These results are independent of the 
choice of bulk potential V{^) for the scalar. Any dependence on V comes in starting at the 
next order in Vq. (2) It may seem strange that on the Planck brane (y = 0) powers of 
track powers of p*, but not so on the TeV brane. The difference is due to the corrections 



to ^ from 6Ai ~ H6Ai, eq. ( p.l3|) , and 6Ni, eq. (p!4|) . They give contributions to of 
order x [6Ai{l),6Ni{l)]. Since ~ /)* + p at lowest order, the correlation of powers 
of 17^ and powers of p* is no longer respected. 

Equipped with these results, we can now specialize to several situations of cosmological 
interest. 



4.1 Inflation 

The earliest epoch of interest is an inflationary one, where the equation of state is a; = — 1 
or w^, = — 1, depending on whether we put the inflaton on the TeV or the Planck brane. 
Let us first suppose the inflaton is on the the TeV brane, so that p* = 0. Then eq. (^^) 
tells us that there is no quadratic correction to the Hubble rate at y = 1, in contrast to 
the situation in pure RS-II (single brane) cosmology. Not only is this true for the leading 
term in the expansion in rrir, but in fact we have checked that it is true order by order for 
the terms 1/m^ and m^. The fact that many terms must seemingly miraculously cancel to 
achieve this suggests that it is true to all orders. 

The analogous statements hold true for an inflaton on the Planck brane: if p = and 
only is nonzero, there is no quadratic correction to at y = 0. Again, we have shown 
this to be exactly true for the first few orders in an expansion in m^. There is such a 
correction to at the TeV brane, but unfortunately this has no impact on the rolling of 
a scalar field on the Planck brane, since its own Hubble rate is unaffected. Thus we cannot 
take advantage of the modification to the Friedmann equation to admit steeper inflationary 
potentials than are normally allowed by the slow-roll conditions, as was explored in refs. 
fs], ^. Only if both p^ and p are nonzero do we get this kind of effect (through the p^p term 



in eq. (O)). But even there, its sign tends to be the wrong one for the desired purpose. 



For example, suppose that some of the inflaton's energy density has been converted energy 



on the TeV brane, with an arbitrary value of uj while = —1. Then (4^) predicts that 
is decreased by the quadratic corrections, regardless of the value of to. 
Even though there might be no presently observable effect, it is still of theoretical 
interest to consider the modification to the TeV brane's Hubble rate due to a Planck brane 
inflaton, to contrast with simpler models of brane cosmology. With p = but p^, ^ and 
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H'^\y=i is corrected by the factor 



SttG f 327rG 



ff^u = ^p.^i + i;^f.j (4.10) 

The correction becomes important at temperatures of order the intermediate scale T ~ 
■sJnirMy ~ 10^'^ GeV. We can compare this to the situation where there is no stabilizing 
scalar field, but a fine-tuning p = —p^./Q'^ is required between the energy densities on the 
two branes, in order to keep the bulk from expanding. In this case, the relation (|l.lD applies 
to the given order at both branes, since for = —1 the difference -ff^(O) — H^{1) = 0{pl) 



1 11]. But this correction only becomes important when T ~ Mp in the unstabilized case. 



since the natural scale for the Planck brane tension is M^. 

The effect in the present case is easy to understand in terms of the usual slow roll 
condition for scalar fields during inflation. If the Hubble rate exceeds the mass of the 
radion, it will start rolling, hence the size of the extra dimension is destabilized, which in 
turn changes the strength of gravity from the 4-D point of view, due to the dependence of 



G on 6o in ( 3.17] ) . The criterion that H < translates to T < y^rrirMp, which is precisely 
the condition we found for the pi correction to to be small. On the other hand, the 
reader should not be misled into thinking that the pi term can be deduced from this effect 
(Hubble expansion giving an m^-dependent shift in b, leading to a shift in G alone). The 
coefficient of pi is numerically different from that which is obtained solely from the effect 
of shifting G. 

We noted that an inflaton on the TeV brane gives no fP correction to measured 
on the TeV brane. A similar story applies for an observer on the Planck brane: he sees 
no p1 correction to on his brane. Again, we find this to be true order by order in an 
expansion about = 0, so it is presumably an exact statement. Moreover we have found 
this to be true for either of the choices of the bulk potential V{^) that we have considered. 
At first sight this may be a surprising result, because it holds even in the limit that the 
infrared brane is removed to y — > cxd. In this limit we might expect physics to coincide with 
that of the RS-II model, where the pi correction to H'^\y=o (eq. (|l.lD ) is well known to 
occur, at least in the case of no bulk scalar field. We will explain the apparent discrepancy 
in a subsection below. 

4.2 Radiation era 



Our results ( [4.7[j4.9| ) predict that there is no effect of order p^/m^ when the equation of 
state on the two branes is that of radiation, to = a;^= = 1/3. The origin of the 1/m^ 
dependence is the large shift in 5bi ( p.l| ) for generic equations of state. In the special case 
of radiation, 6bi = since the radion couples to the trace of the stress energy tensor, which 
vanishes in this case. We have therefore done a separate treatment when a; = = 1/3 to 
find the leading effect. We obtain 

HX=i = ^ ((p. + P) + ^,i^'p* - P)P*) (4.11) 
2, _ SttG ( 2ttG 2 



H%=o = ^{ip*+p)- ^j^i^ P* - P)P] (4-12) 
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and the equations for dH/dr have vanishing corrections at this order in p and p-t. It is 
curious that, just Uke in the previously discussed case of inflation, the Hubble rate on a 
given brane receives an order p^ correction only if the energy density on the other brane is 
nonvanishing. 

We can make the preceding results much stronger in the case where the energy density 
vanishes on one of the two branes. If = 0, then eq. (3.13) implies that the second order 
correction to the Hubble rate on the TeV brane vanishes identically, with no approxima- 
tions, provided that 



p4Aofei) 



~2Ao(yi) 



-2A(,fe) 



4^0(1/2) 



-2Ao(?;3) 



„4Ao(?/i) 



yi 



y2 



-2Ao(s/3), 



dys dy2 



dyi 







(4.13) 



In the inverse case where p = and cj* = 1/3, the second order correction to evaluated 
on the Planck brane vanishes if 



+ e 



-2Ao(j/i) 



e 
ryi 



-2Ao(j/2)^y2 
,4Ao{j/2) 



-2Ao(?/2) 



dy^. 



(4.14) 



J^' e^^^^y^^ (^j^J e-^^^^y^Uy^ - ^ e-^^'^^y^^ dy^ dy2 



dyi 



Writing 



Ul 

dv 



yi 



yi 



Mo(y2) 



=4Ao(s/2) 







y2 



y2 



-2Aofe) 



'•''dy3 dy- 



-2Aofe)^, 



2/3 



-2Ao(y3) 



'•"'dyz dy2 



-2Ao(2/i) 



dyi 



(4.15) 



where ui and U2 are to be used in the flrst and second case respectively, we can use 
integration by parts (J uidv = UiV — J vdui) to show that these integrals do indeed vanish 
identically. 

Let us now focus on physics on the TeV brane {y = 1). By the time of nucleosynthesis 
no more than 10% of the radiation energy density can be on the Planck brane, so one 
expects that p^ ^ p at TeV temperatures as well. Therefore we can ignore the Q.'^p^, term, 
and the fractional correction to is of order 



5H^ 



1 



( P* 
VTeV^ 



(4.16) 



It is interesting that the order of magnitude is nearly correct for this to be relevant at 
the time of the electroweak phase transition. (We could change TeV 100 GeV by 
making k a few orders of magnitude smaller than Mp.) However the minus sign shows that 
the Hubble rate is suppressed relative to standard cosmology. This is disappointing for 
electroweak baryogenesis, since one would have liked to increase make sphalerons go 
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out of equilibrium at the electroweak phase transition, a possibihty that has been recently 
analyzed in ref. ||2^. Of course, it may happen that the sign turns out to be the right one 
in some other model of brane cosmology. 

4.3 Relation to RS-I model without stabilization and to RS-II 

As a check on our results, one would like verify that they are consistent with related 
cosmological solutions. For example, when there is no scalar field and the equation of state 
is inflationary, Kaloper |27] found the exact solution^ 

a{w) = e'^"" - sinhku! (4.17) 

where p^^ = p^^/T, the ratio of the excess energy density to the tension on the Planck brane, 
and the Hubble rate is given by 

= 2k'^{{l + p^f - I) (4.18) 

As in the static case, one can decide whether or not to compactify the extra dimension 
by inserting a second brane at some position wi. Without compactification, the extra 
dimension is still effectively cut off by the presence of a horizon at = ^ ln(l + J-). If 
we compactify explicitly by inserting a second brane, its tension (— T) and excess energy 
density p are completely fixed by the value of a'{wi), leading to the constraint 

J e '^^^ — smh KWi 



Expanding in p^, to leading order this gives the well-known constraint |11| p = —p^/fl'^ 
needed for a static bulk in the absence of stabilization (recall that the Q = e~^'^^ is the 
warp factor). However, the exact tuning is a power series in p^,, and our second order 
formalism correctly reproduces the 0{p1) term, as we describe in more detail in appendix 
C. Of course, our perturbative treatment also correctly reproduces the solution a{w) and 
ao{t) expanded in powers of p*. 

The previous check was specifically for the equation of state = —1. The same can be 



done for an arbitrary equation of state, for which the (p^l) solution generalizes to |25, 29| 

ds'^ = -n^{t, w)dt^ + a^{t, w)dx'^ + dnP'] 
a{t,w) = ao{t){e~'^'" - p^.{t)sinhkw); 

n{t,w) = e~'"" + p^{t){2 + 3uj^) smhkw, (4.20) 



where oo(t) is determined by the modified Friedmann equation (4.18) with H = Again, 
our perturbative formalism agrees with the expansion of the exact solution at order pi . It 
is interesting to note in passing that for > —1, the metric function a{t,y) vanishes at a 



*This is related to the coordinate system we have been using by w = by, where now b is taken to be 
constant by fiat. 
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position between the brane at y = and the horizon where n(t, y) = 0. Nevertheless, all 
components of the 5D Riemann tensor are well behaved throughout the bulk, and neither 
the vanishing of a{t, y) nor n{t, y) leads to any curvature singularities. 

Now we come to a slightly more subtle issue: what happens if we try to smoothly 
remove the TeV brane from the problem, to recover the known single-brane results? We 
have noted that our results for do not reduce to those of the RS-II model, even in the 
case where the TeV brane is taken to infinity. The presence of the bulk scalar is crucial 
for understanding the difference between the two situations, because it modifies the bulk 
stress-energy component T55 in response to matter on the branes. Even though we are 
working in a gauge where 5<^i = 0, T55 = 555(V^(^) + A) — is nevertheless modified by 
the shift in g^^ = b{t,y)'^. It is the G55 = k?T^5 Einstein equation which implies the exotic 
form of the Friedmann equation ~ (p* +T)^ when there is no matter in the bulk. This 
argument no longer holds when T55 is present. Therefore it is not so surprising that the pi 
corrections no longer coincide with those of the RS-II model, when we have just a single 
brane but also a scalar field. 

Mathematically, we cannot recover the RS-II model simply by taking the limit vq — > 0, 
which corresponds to removing the scalar field from the problem. If we attempt to do 
this, our perturbation expansion breaks down because 6bi, eq. ( |3.1[ ), diverges — unless the 
matter on the branes is tuned so that the term in brackets in ( |3.1| ) vanishes. This tuning 
is the condition p = —p^/Q,'^ discussed above. A special case is that of radiation, where no 
such breakdown of perturbation theory occurs. 

4.4 Dark Radiation 



In the present work we have neglected the presence of dark radiation |28, 3^, 31, 32], which 
appears in the Friedmann equation in the normal way, 

i/^ = 2.'^((p. + lf-l + ^); (421) 

in particular, there is no quadratic correction /ao{t)^ even in the unstabilized case. The 
AdS/CFT interpretation is that the dark radiation is simply the thermalized degrees of 



freedom of the conformal field theory |17], but from the 5D viewpoint, the dark radiation 
is associated with the presence of a black-hole-like singularity that appears in the bulk at 
a finite distance from the brane. This can be most clearly seen in the coordinate system of 
ref . (see also ) , where the brane is moving through a static bulk whose geometry is 
that of the 5D AdS-Schwarzschild metric. 

Again using the results of [^8|, the exact solution takes the form ([4.20| ), but now 



a\t,w) = al{t) (e-'=--p,(t) 



2 



sinh kw ] H ^-^tj sinh^ kw; 

aoity 



n{t,w) = ^ = ^ ( (e->^- -Mt) sinh kwY (4.22) 
ao{t) a{t,y) \\ J 

+3/9*(t)(l + Ld*) (e~^^ — p^:(t) sinhkw^ sinh kw 3^ sinh^ kw] 

V J a^{t) J 
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From this exact solution we see that the metric functions do not factorize in the way which 
we assumed in our original ansatz ( |2.3| ) . A more general form would therefore be needed to 
study the interactions, if any, between corrections and dark radiation in the stabilized 
model. 



5. Summary and conclusions 

In this paper we have focused on the corrections to the Friedmann equation which arise in 
the possibly more realistic 5D brane cosmologies where the bulk stress energy is not trivial, 
but has some inhomogeneity in the extra dimension due to a scalar field. The presence of 
the scalar was motivated by the need to stabilize the size of the extra dimension when the 
extra dimension is compactified by the presence of a second (TeV) brane, but it could also 
be present on more general grounds. 

In a perturbative expansion in powers of energy densities on the branes, the Hubble 
parameter H generally should receive corrections to all orders in p and /)* rather than 
terminating at order as happens in the unstabilized solutions. We found the rather 
surprising result that the correction vanishes in the special cases where the energy 
density is confined to the same brane as that where H is measured and the equation of 
state is either that of inflation or radiation. If either of these conditions are not fulfilled, 
the p^ correction does not vanish, but its sign and magnitude can be different from the case 
of a single brane with no bulk scalar field. In these cases, the correction has a coefficient of 
order (TeV)~^m~^, where is the mass of the radion, which is expected to be somewhat 
lighter than the TeV scale. 

We noted two implications of this result. First, it is impossible to obtain order p^ 
corrections to during inflation on a given brane unless the inflaton is on the opposite 
brane. This means the corrections have no effect on the dynamics of the inflaton, which 
is disappointing from the point of view of potentially observing the effects of branes on 
cosmology. Second, it is equally true during radiation domination that p^ effects can appear 
for a given brane only if the energy density on the other brane is nonzero. If some fraction of 
the radiation of the universe is on the Planck brane, there could be significant deviations 
from standard cosmology during the electroweak epoch. Unfortunately, the sign of the 
corrections is such as to decrease the Hubble rate. It would have been more interesting for 
the purposes of electroweak baryogenesis to obtain the other sign. 

On a more optimistic note, we should reiterate that the exact cosmological behavior 
cannot be predicted from the p^ terms when they become as important as important as the 
linear in p terms, since the perturbative expansion is breaking down. It would therefore be 
interesting to find exact cosmological solutions in the presence of the stabilizing field. 

We thank Hassan Firouzjahi for enlightening discussions, and Geraldine Servant for 
helpful comments on the manuscript. 
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A. 1st order functions appearing in 2nd order Einstein equations 



These are the functions appearing in eqs. which give the solutions for the 

second order perturbations. 



55 



8hi 
^0 



- - (^*'^- 4) 4 - ^*^') 



+^ y \^^\ - ci>;,^ j + (^ci>'o^ - Jc^l) j (A.3) 



.2, 

+y UTi - + 4^1)^0 + (-2-^ - 45iVi)^o + -^^-'I'J 5<i>i (A.4) 

Oo Oo Oo 

+5'^[ + + 4^1 + ^-^') + 25<I>'/^ 

V "0 J bo 



B. Complete second order corrections 

In the previous sections, we presented simphfied expression for and H. Here we present 
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the complete expressions for the second order in p corrections to 0{vq evaluated on the 
branes. 

+ (1 - 3a;,)(4(4 + - J72(7 + ?>uj^))pi 

+21^2 |jl2^(i _ 3^)(i _ 3c^^) + 6(1 + w,)(l - 3w,) - 2(1 - 3a;,) - 4(1 - 3w)) 
+2(1-3w) + 6(1+l<j)(1-3l<;)}/9p*] (B.l) 

+0^(1 - 3w)(4172(4 + 3a;) - (7 + 3a;))/92 

+21^2 {(1 _ 3t^)(i _ 3cj^) + 6(1 + a;)(l - 3a;) - 2(1 - 3a;) - 4(1 - 3a;,) 
+f^2(2(l - 3a;,) + 6(1 + aj,)(l - 3a;,))} pp,] (B.2) 

+ (1 - 3a;,)(l + a;,)(4(l - f]^) + 9(1 + ^^))pl 

{(1 - 3a;)(2(l + a;,) + 6(1 + ujf + 2(1 + a;) + 3(1 + aj)(l + w,) 
+02(-2(l + a;)(l - 3a;,) + 2(1 + a;,)(l - 3a;,) + 3(1 + w,)(l - 3w,)(l + lo) 
-4(1 - 3w)(l + w,) + 6(1 + w,)2(l - 3w,))} pp^] (B.3) 

-17^(1 - 3w)(l + w)(4(l - 17^) + 91^2^1 + uj))p^ 

+1^2 _ 3t^^)l72(2(i + ^^) + 6(1 + LO^f + 2(1 + a;) + 3(1 + a;)(l + w,) 

-2(1 + a;,)(l - 3w) + 2(1 + w)(l - 3w) + 3(1 + w)(l - 3w)(l + w,) 

-4(1 - 3a;,)(l + a;) + 6(1 + ujf{l - 3a;)} pp^] (B.4) 

C. Absence of a scalar field and infinite extra dimension 

In order to trust the perturbative approach which was used to derive the results presented 
in this paper, we must convince ourselves that the same approach allows us to reproduce 
the well known results of RS-I and RS-II cosmology in the absence of a scalar field. In this 
section, we will show that this is indeed the case. 

In the absence of a stabilization mechanism, the configuration with two parallel branes 
lying at fixed positions along the extra dimension is generally unstable. In order to look 
for static solutions, we must impose h = constant from the start, which means that any 
derivatives of h are absent from the Einstein equations ( |2.6| ). 

Working to zeroth order in /?, it is straightforward to show that: 

v. = -v. = % 

Aq = kboy, (C.l) 
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so that we recover the expected fine-tuning between the branes' tensions. 

The (00) and (00-ii) combinations of the Einstein equations, hnearized in p, can be 



solved in the same manner as in [15|. However, the absence of a scalar field means that 

(5bi 
bo 



the (55) equation no longer allows one to solve for Rather, plugging the solutions for 



^i,Ti, {^^^ and (^f^^ iiito the (55) equation leads to the following constraint: 

5^^p-3p) +n^(/o* -3p*) = 0. (C.2) 

The only way to satisfy this constraint without running into inconsistencies at higher order 
in p is to set 

p=—Vl~'^p^; uj^ = uj = —I, (C.3) 

where the second equation is found by demanding that ^ (/fr) constant in time, as 
dictated by the first equation. 

In short, setting ^ and 5h to zero has made our system of equations over-determined: 
the (00), (00-ii) and (55) Einstein equations constitute a system of three equation for two 
unknowns. The way out of this problem consists of turning one of our parameters {p and 
/9*) into a variable to be expanded in powers of the other. We will therefore choose to write 
p = Pi + P2 + ■ ■ ■■, where the subscripts indicate the order in powers of , and pi is given 
by Q. 

We can then repeat the same steps (solving the first two Einstein equations and plug- 
ging the results in the (55) equation) for the equations at second order in to find: 

-^2(1-172) 

These results for the relation between p and p^ agree with the exact results presented in 



1 27], when the latter are expanded as a powers series in p^. 

Imposing these constraints we find the Hubble rate on the Planck brane to be 



H\=o = ^ii-n')p. 



+ 0{pt) (C.5) 
with SirG given by ( 3.17]) . Using (CA) and the fact that — > as 5o — > oo, we see that 



we do indeed recover the expected RS-II behaviour (|l.lD when the second brane is taken 
to infinity. We emphasize however that this is only true insofar as pi and p2 are written 
as in ( |C.3| ) and ( p. 4]) . If we had instead decided to take the TeV brane out of the picture 
by setting p = at the start, we would not have gotten this result. Specifically, there 
would be a missing factor of 1/2 in the second order correction, and more importantly, the 
constraint equations ( p.3| , p.4| ) would demand that = 0. 

Suppose now that we start out with a single brane with an extra dimension of infinite 
size. The difference between this approach and the one described above lies in the fact that 
setting p = also sets the boundary conditions for ^ and T at y = 1 to be ^'(y = 1) = 
T(y = 1) = 0. Here however, since the second brane is explicitly absent from the setup, 
we only have one set of boundary conditions on ^ and T, i.e. these two variables are only 
fixed on the single brane, located at y = 0. 
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We start by redefining the extra dimension's coordinate y as 

y = boy, (C.6) 

so that ij goes from to 60 as y goes from to 1 (and 60 can be chosen to be infinite). We 
also choose a gauge where the fluctuations in h vanish. To order , the equations tell us 
that 

^0 = ky 



The order Einstein equations will again look like those of |15|. Even though we only 
have one set of boundary conditions for ^1 and Ti, we can still solve the first two Einstein 
equations, leaving ^i{y = bo) and Ti(y = 60) unspecified. If we then plug the results in 
the (55) equation, we find the following constraint: 

(4^1 + Ti) \^=to = ^'^ im + Ti) 1^=0 (C.8) 
Taking this into account, we can write the Friedmann equations as: 



Hi = Q^^^^2^ + 317V - 3 - 3uj,)p, + 2in^^i\^=bo) (C.IO) 

It would appear that we have complete freedom in choosing a value for ^i{y = bo), but 
there are a couple of points to consider. 

Firstly, we expect our results to reproduce the standard Friedmann equations at 0{pl). 
One can easily see that with the following choice: 

*i(y = M = ^p* (c.ii) 

D 

this requirement is satisfied. Similarly, recovering the correct 0{p1) term would require 
setting ^'2 to the appropriate value at y = 6o- 

Secondly, we know from [^] that if the tension on the brane is larger than the static 
solution's value, there will be a horizon in the bulk at the (finite) value of y for which 
^{y) = 0- Given the form we have chosen for a, it is not immediately obvious how we 
can reproduce this behaviour. However, if we expand e~^^^ as 1 — 5Ai, then in the limit 
where we neglect any higher order in p contributions, a will vanish when 5Ai = 1. (Notice 
that this is also the point at which our perturbative approach breaks down completely). 
Choosing ^i{y = bo) in the manner mentioned above, and recalling that 

Mi(y) = f\i{y)dy (C.12) 
Jo 

we find that 5Ai{y) = 1 at the position 
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which corresponds precisely to what we would find using the exact solutions of |27]. 

One more subtlety remains. If we solve for 6Ni, we can see that there will be a point 
in the bulk where n vanishes. (In the same manner as o vanishes, i.e., where 6Ni{y) = 1). 
In general, this point and the one where a vanishes do not coincide. Indeed, it can be 
shown that 6Ni = 1 at the position 

1 / 12k \ 
yn=o = ^11^1 + ^ + 3(1 + - 1) J ■ (C-14) 

Intuitively, one might expect that a = corresponds to a geometrical singularity while 
n = corresponds to a horizon. We would then demand that y„=o < ya=o in order 
to avoid the appearance of a naked singularity in the bulk. One can readily see that this 
would impose < —1. However, as we stated in the main text, examining the behaviour of 
the Riemann tensor at these special values of y shows no evidence of curvature singularities 
in the bulk. 

To summarize, we have shown in this section that our perturbative approach allows us 

to: 

a) reproduce the expected RS-II result as a limit of unstabilized RS-I, as long as p is 
treated as a dependent rather than a free parameter; 

b) reproduce the expected RS-II results starting from a single brane setup, as long as 
we fix the values of ^ and T at y = 60 appropriately; 

c) calculate the points in the bulk at which a and n vanish and find a result which 
matches the exact one. 
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